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Abstract 

We construct extensions of the standard model based on the hypothesis that 
the Higgs bosons also exhibit a family structure, and that the flavor weak 
eigenstates in the three families are distinguished by a discrete Zq chiral sym- 
metry that is spontaneously broken by the Higgs sector. We study in detail at 
the tree level models with three Higgs doublets, and with six Higgs doublets 
comprising two weakly coupled sets of three. In a leading approximation of 
S3 cyclic permutation symmetry the three Higgs model gives a "democratic" 
mass matrix of rank one, while the six Higgs model gives either a rank one 
mass matrix, or in the case when it spontaneously violates CP, a rank two 
mass matrix corresponding to nonzero second family masses. In both models, 
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the CKM matrix is exactly unity in leading approximation. Allowing small 
explicit violations of cyclic permutation symmetry generates small first family 
masses in the six Higgs model, and first and second family masses in the three 
Higgs model, and gives a non-trivial CKM matrix in which the mixings of the 
first and second family quarks are naturally larger than mixings involving 
the third family. Complete numerical fits are given for both models, flavor 
changing neutral current constraints are discussed in detail, and the issues of 
unification of couplings and neutrino masses are addressed. On a technical 
level, our analysis uses the theory of circulant and retrocirculant matrices, 
the relevant parts of which are reviewed. 

I. INTRODUCTION 

It has long been recognized that the hierarchical structures of the family mass spectra, 
with their large third family masses, and of the CKM mixing matrix, with its suppressed 
third family mixings, may have a common dynamical origin. In particular, several authors 
[1] have stressed that the observed pattern seems to be close to the "rank-one" limit, in which 
the mass matrices have the "democratic" form of a matrix with all matrix elements equal 
to unity, which has one eigenvalue 3 and two eigenvalues 0; when both up and down quark 
mass matrices have this form, they are diagonalized by the same unitary transformation 
and the CKM matrix is unity. A generalization of the democratic form, that is closely 
related to the models developed below, is the suggestion of Harrison and Scott [2] that the 
Hermitian square of the mass matrix should have the form of a circulant matrix. Because the 
underlying dynamical basis for these choices has not been apparent, it has not been possible 
to systematically extend them to renormalizable field theory models that incorporate, and 
relate, the observed mass and mixing hierarchies. 

We present in this paper models for the quark mass and flavor mixing matrices, based on 
the underlying dynamical assumption that the three flavor weak eigenstates are distinguished 
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by different eigenvalues of a discrete chiral quantum number. The idea that a discrete 
chiral quantum number may underhe family structure was introduced originally by Harari 
and Seiberg [3], and was developed recently by the author [4] in a modified form that we 
follow here. Also of relevance is the remark of Weinberg [5] that an unbroken discrete chiral 
quantum number suffices to enforce the masslessness of fermionic states. Extending the 
general framework of this earlier work, we postulate that all complex fields carry a discrete 
chiral family quantum number. Since the Higgs scalars in the standard model are complex, 
we introduce one or two triplets of Higgs doublets that carry quantum numbers, and 
that are coupled to the fermions by Yukawa couplings constructed so that the Lagrangian 
is exactly Z^ invariant. Spontaneous symmetry breaking, in which the neutral members of 
the three or six Higgs doublets acquire vacuum expectations, then gives the fermion mass 
matrices that form the basis for our detailed analysis. 

In addition to postulating that the Lagrangian has an exact discrete chiral symmetry 
that is spontaneously broken, we also postulate that there is an Sj, cyclic symmetry under 
cyclic permutation of the flavor eigenstates that is explicitly but weakly broken by the 
Yukawa couplings and the Higgs self-couplings in the Lagrangian. This assumption permits 
the analysis of our models by developing them in a perturbation expansion in powers of 
the ^3 cyclic symmetry breaking, leading, as we shall see, to qualitative features of the 
mass and mixing hierarchies that accord with observation. An interplay of spontaneously 
broken symmetries with weakly explicitly broken symmetries has played a useful role in 
particle phenomenology in the past, most notably in understanding the consequences of 
chiral symmetry in quantum chromodynamics. Our analysis suggests that such an interplay, 
in the context of electroweak symmetry breaking, may also provide a basis for understanding 
features of the mass and mixing hierarchy. 

This paper is organized as follows. In Sec. II we elaborate on the form of and motivation 
for our basic assumptions of an exact discrete chiral symmetry and an approximate S-i, cyclic 
permutation symmetry. In Sec. Ill we write down the Lagrangians for two extensions of the 
standard model that incorporate these assumptions, the first based on a single three family 
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set of Higgs doublets, and the second based on including an additional weakly coupled 
three family set of Higgs doublets. In Sec. IV we review the theory of circulant and 
retro circulant matrices, in the framework of the 3x3 matrices that are needed for the 
subsequent analysis. In Sec. V we discuss the extrema of the Higgs potentials in the three 
and six doublet models, in the limit of exact S3 cyclic symmetry. We work out the spectra of 
physical Higgs particles, and show that for a wide range of parameters, the six doublet model 
leads to spontaneous violation of CP. In a related Appendix we give the formulas needed 
for numerical minimization of the Higgs potentials by the conjugate gradient method. In 
Sec. VI we use the extrema determined in Sec. V to calculate the tree approximation mass 
matrices. We show that in the limit of exact cyclic permutation symmetry, the mass matrices 
are retrocirculants, corresponding to the rank one "democratic" form in the three doublet 
model and to a rank two generalization in the six doublet model when CP is spontaneously 
violated. Also, in the limit of exact cyclic permutation symmetry, we characterize the Higgs 
decay modes, and show that the CKM matrix is exactly unity and that strangeness changing 
neutral currents exactly vanish. In Sec. VII we formulate a perturbative expansion around 
the zeroth order approximation of exact Ss cyclic permutation symmetry, and show that the 
mixing matrix for the first and second families is zeroth order in the perturbation, whereas 
the mixings involving the third family are first order in the perturbation. In Sec. VIII we 
derive formulas for the contributions from Higgs exchange to the Kl — Kg mass difference, 
which is the process most sensitive to strangeness changing neutral current effects. In Sec. 
IX we describe the procedure used for making overall fits of our model, including small 
violations of cyclic permutation symmetry, to the data, give sample numerical results, and 
draw some conclusions from these. In Sec. X we summarize experimental signatures for our 
model, comment on its extension to neutrino masses and mixings, discuss the prospects for 
coupling constant unification, and give some directions for future investigations. 



4 



II. BASIC ASSUMPTIONS: AN EXACT DISCRETE CHIRAL SYMMETRY 
AND AN APPROXIMATE ^3 CYCLIC SYMMETRY 

In formulating our basic assumptions, we shall follow a procedure that has worked well in 
the past as an heuristic tool in particle physics. This is to abstract symmetry or partial sym- 
metry assumptions from specific simplified field theory models, then to discard the models, 
but to retain the symmetry assumptions deduced from them as the basis for phenomeno- 
logical calculations. Examples where this has been a productive method in the past include 
(1) the CVC and PCAC symmetries of the strong interactions, the algebra of currents, and 
the calculational methods based on these, and (2) the approximate SU (3) flavor symmetry 
of the strong interactions. These postulates, which had a somewhat ad hoc character at 
the time when they were first introduced, helped pave the way for the formulation of the 
standard model, into which they were incorporated in a natural way and thereby ultimately 
justified. 

Our aim in this paper is to apply a similar method to the problems of family structure 
and mass and mixing matrices, which to date have been among the most vexing puzzles 
of the standard model. As an heuristic field theoretic model, we shall adopt a simplified 
composite model in which all matter particles (quarks, leptons, and Higgs fields - every- 
thing other than the gauge fields) are composites of a single fermion field x- As observed by 
Harari and Seiberg [3] and Weinberg [5] , in a gauge theory for x the instanton determinant 
that breaks global U{1) invariance leaves unbroken a discrete Z'2.k chiral subgroup, with 
K determined by the index of the representation of the gauge group under which x trans- 
forms. Harari and Seiberg propose, moreover, that this naturally occurring discrete chiral 
subgroup provides the quantum number that distinguishes between the various families. 
Since it is now clear that there are exactly three light families, we shall assume henceforth 
in applying this idea that i^T = 3, so that we start from the assumption that the funda- 
mental Lagrangian, as augmented by the instanton-induced potential, is invariant under the 
simultaneous transformations 
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XL ^ XL exp(27ri/6) , ^ Xi? exp(-27ri/6) (la) 

of the fundamental fermion fields x- The fields in the low energy effective Lagrangian are 
in general nonlinear functionals of the fundamental fields. Fermionic effective fields must be 
odd monomials in the fundamental fields, and so can come in three varieties ijjn with the 
discrete chiral transformation law 

V'nL VnLexp((2n + l)27ri/6) , VnR '0nflexp(-(2n + l)27ri/6) , n = 1, 2, 3 , (lb) 

while complex bosonic effective fields must be even monomials in the fundamental fermion 
fields, and so can also come in three varieties 0„ with the discrete chiral transformation law 

0n — 0nexp(2n27ri/6) , n = l,2,3 . (Ic) 

Introducing the cube roots of unity u and ZJ, 

1/3 1/3 
uj = exp(27ri/3) = "2 + ^ = exp(-27ri/3) = -3 ~ ' C^^) 

that obey the relations 

Zu = cu* =cu'^ , l+a; + aJ = , (2b) 
the transformation laws of Eqs. (la-c) take the form 

XL Xluj^ , Xr XRM^^ , 
IpnL i^nLi^''^^ , l/jnR ^ IpnR^^"'^^ , n = 1, 2, 3, 

0„^</.,(^^ n = 1,2,3 . (3) 

Gauge fields are real fields, and since the phase in Eq. (Ic) never takes the value —1 for 
any n, the gauge fields in a Zq model necessarily come in only one variety transforming with 
phase unity under discrete chiral transformations. Thus the minimal Zq invariant extension 
of the standard model consists of a triplicated set of fermions, and a triplicated set of Higgs 
doublets, obeying the transformation laws of Eqs. (lb) and (Ic) respectively, together with 
the usual gauge bosons, with the Lagrangian constructed to be Zq invariant. 
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As we shall see in Sec. Ill below, the assumption of an unbroken discrete chiral symmetry 
still leaves many parameters in the Lagrangian, and it is desirable to look for a further exact 
or approximate symmetry to impose. The natural candidate is 5*3 cyclic permutation sym- 
metry, under simultaneous cyclic permutation of the n — 1,2,3 discrete chiral components 
of the fermion and Higgs boson fields. If the discrete chiral components were physically 
identical, one would expect this 5*3 cyclic symmetry to be exact. However, in the compos- 
ite picture from which we are abstracting our model, the discrete chiral components differ 
physically by the addition of fermion- antifermion pairs coupled as Lorentz scalars, and so 
the internal wave functions of the discrete chiral components are different. Thus the best 
we might hope for is an approximate, weakly broken, 5*3 cyclic permutation symmetry, and 
this will be assumed as the second ingredient of our model. 

By abstracting our two fundamental assumptions from a schematic composite model, 
we gain some assurance that they are consistent with each other and at least physically 
plausible. However, we do not attach great significance to the particular model from which 
they were inferred; it is entirely possible that the same assumptions can emerge from other 
dynamical frameworks. We shall henceforth avoid further discussion of underlying models, 
and focus on exploring the consequences of our assumptions within the standard framework 
of low energy renormalizable effective action phenomenology. 



III. DISCRETE CHIRAL INVARIANT EXTENSIONS OF THE STANDARD 

MODEL 

We proceed now to write down discrete chiral invariant extensions of the Lagrangian 
density for the standard model, following the notation of the text of Mohapatra [6]. In the 
following, each quark or lepton field is implicitly a column vector formed from the three 
discrete chiral components obeying the transformation laws of Eq. (3), with the n — 1 index 
at the top of the column vector and the n = 3 index at the bottom. For the Higgs scalar 
fields, the discrete chiral subscript n will be indicated explicitly. We shall be interested 
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in two models, the first containing a single discrete chiral triplet of Higgs doublets 0, the 
second containing two discrete chiral triplets of Higgs doublets, denoted respectively by 
and T]. We shall write all formulas for the case of the six Higgs doublet model; the simpler 
three doublet model is obtained by setting all fields 77 to zero. 

The total Lagrangian density £, consists of kinetic terms for the gauge, Higgs, and 
fermionic fields, together with Yukawa couplings of the Higgs fields to the fermions and 
the Higgs self-interaction potential. Writing 

^ — ^gauge kinetic ~l~ ^Higgs kinetic ~l~ ^fermion kinetic 
Yukawa ~l~ -^Higgs potential ) 

(4a) 

the gauge kinetic terms have the usual form 

■^gauge kinetic "^^^lu ' ^^^lu ~^B^yB^y , 

W^,^d^W,-d,W^ + gW^xW, , B^,^d^B,-d,B^ , (4b) 

and so also do the fermion kinetic terms (with Ql and -0^ respectively the left-handed quark 
and lepton doublets, and r the weak isospin Pauli matrices that act on them), 

Aermion kinetic = ~Q lIA^IJ- ~ '^^ ' ^IJ- ^B^)Ql 

-V^l7.(5m-|^-W>;.+ |^Wl 

-UR-i^id^ - '^Bi^)ur - dR'y^{df, + y 5^)c?fl . (4c) 

The Higgs kinetic energy is simply a sum over kinetic terms of the standard form for the 
discrete chiral components of the scalars and 77 (each of which is, as usual, a weak isospin 
doublet), 

-^Higgs kinetic = — X] \(^tJ.4>n — ■ W ^(pn 

n=l,2,3 ^ ^ 

- E \d,Vn-'-^T-W,r^r,-'-^B,r^^\' . (4d) 

n=l,2,3 ^ ^ 



It is only in the Yukawa couplings and the Higgs potential that invariance under discrete 
chiral transformations plays a nontrivial role. Letting 0„ and fin denote the CP conjugates 
of the Higgs fields, 



0, = (CP)- VnCP = ^T2C 



(5a) 



the Yukawa Lagrangian takes the form 



■'-'Yukawa — 

Ql^'^cIr + Ql^'^ur + ^i$"eR + + adjoint, 



(5b) 



where , / = (i, M,e,z^isa3x3 matrix acting on the discrete chiral column vector structure, 
and where we have allowed for the possibility of nonzero Dirac neutrino masses by including 
a right-handed neutrino. The matrices must be constructed so that Eq. (5b) is invariant 
under simultaneous discrete chiral transformations of the fermion and Higgs fields. Referring 
to Eq. (3), it is easy to see that this dictates the structure. 



¥ = giiP^i + + P/303) + 4(/^^i + + /^%), f = d,e , 



(6a) 



pf 



pf _ 



1+A 








/ 

■^21 







with the 3x3 matrices P^^. given, for all flavors f — u,d,e,v and for ^ = 0, t;, by 

/ 1 + /3/12 \ 






1 + (^{22 





1 + f^u 

To uniquely specify the Yukawa couplings gf^, we require that the parameters sum to 
zero. 



pf _ 





V 










^23 



(6b) 



/ 



E/5inn = . (6C) 
mn 

When there is exact S3 cychc permutation symmetry the /3's all vanish, and thus the case 
of approximate 5*3 cyclic symmetry is parameterized by /3's that are all small compared to 
unity. In a CP conserving theory all of the coupling constants g^^^ and all of the /3's are 
real; when CP conservation is not imposed, these parameters can be complex. 

We turn finally to the Higgs potential, which we separate into four terms as follows, 

^Higgs potential = + + 7^) + V2(0, 7]) , (7a) 

with (for ^ = (p,r)) 

n=l 

V^n = >^^niCiCn " V^J^ - Hl^nClCnd+lCn+l " l^2^n\CiCn+lf " a^n^e exp(iV'^„)^j,^„+i^j,^„_l , (7b) 

where the coefficients in Eq. (7b) are real (by hermiticity) and where the parameter ip^n is 
zero (modulo tt) when CP conservation is imposed. For the potential terms that couple the 
4> and rj Higgs fields, we have in the CP conserving case 

3 

m,n=l 

+C3mnRe(l>ln(l^m+lVliVn-l + (^47^71 Rer^m^m+l^Ji^n-l 

+C5mn^(i<PlnVm+lVi4>n-l + C'emnReT^^^m+l^l^n-l) , 

^2(0, ?7) = Y^nRe(l)iVn 
n 

3 

m,n=l 

+CQmn^er}l^Vm+l<PlVn-l + CiOmn'Rer]l^'nm+lVi4>n-l 

+CiimnRe(l>lnVm+l(l>liVn-l + C'l2mnRe77j„0m+l?7n0n-l) , (7c) 

with all constants real (again by hermiticity). The terms Vi are those invariant under 

independent rephasings (pn Gxp{i9if,)(f)n and rjn exp{i9r^)rin of the two Higgs discrete 
chiral triplets, while the terms V2 are only invariant under this phase transformation when 
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restricted so that — 9jj. When CP is not conserved, an independent phase can be inserted 
inside each real part Re in the above expressions, in analogy with the construction of the 
final term of Eq. (7b). When there is S3 cyclic permutation symmetry, the constants with a 
single discrete chiral subscript n are independent of that subscript, while the constants with 
a double subscript mn obey the cyclic condition Cimn = Ctm+in+i, ^ = 1, 12 . 

This rather complicated Higgs potential completes the specification of our model, the 
tree approximation to which will be analyzed in detail in the sections that follow. 



IV. RETROCIRCULANT AND CIRCULANT MATRICES 



Before proceeding further, we pause to review the theory of circulant and retro circulant 
matrices in the 3x3 case relevant for what follows. For a compact summary of general 
results see Marcus [7] and Hamburger and Grimshaw [7], and for a detailed exposition see 
Davis [8]. A matrix 

/ a b c\ 



Circ^(a, b, c) = 



(8a) 



Circ^(a, b, c) 



(8b) 



cab 
\b c a J 

is called a circulant, while a matrix 

/a b c\ 

b c a 
\c a b J 

is called a reverse circulant or retrocirculant. [Clearly, a retrocirculant is always a symmetric 
matrix, and so Circ<_(a, 6, c) — Circ<_(a, 6, c)"^, and Circ^(a, 6, c)^ = Circ^(a, 6, c)*.] Two 
properties of these matrices are used in what follows. The first is that the Hermitian square 
of a retrocirculant is a circulant, 

Circ^(a, b, c)Circ^(a, b, 0)"^ = Circ^(|a|^ + + |c|^, ab* + be* + ca*, ac* + ba* + cb*) , 
Circ^(a, b, c)^Circ^(a, b, c) = Circ_^(|a|^ + + |c|^, a*b + b*c + c*a, a*c + b*a + c*b) . (9) 
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The second is that any retrocirculant with arbitrary complex a,b,c is diagonahzed by trans- 
formation from the left and right by unitary matrices Ul, Ur — Ul, that are independent 
of the values of a, 6, c. Explicitly, setting 





1 


/I 








1 


CO 


w 






\ ^ 


1 

1 


1 , 






/I 


CO 


W\ 




1 

73 


1 


UJ 


UJ 






u 


1 


1/ 






/I 


1 


1\ 




1 

73 


UJ 


UJ 


1 






u 


UJ 


1/ 



(10a) 



a simple calculation shows that 



C/LCirc^(a,6, c)C/i = 












V 



(10b) 



a -\- ujh -\- ujc 

a + b + c) 

An elementary corollary of these statements is that any Hermitian circulant matrix 
is diagonahzed by the unitary transformation UlH^U]^ using the unitary matrix Ul ol 
Eq. (10b). 

The relevance of these results to what follows is that in the limit of 5*3 cyclic permutation 
symmetry, we shall find that the fermion mass matrices in both the three and six doublet 
models arc rctrocirculants, and so are diagonahzed by the universal bi-unitary transfor- 
mation of Eq. (10b). By Eq. (9), the Hermitian squares of the fermion mass matrices in 
the approximation of cyclic permutation symmetry are therefore circulants, as suggested by 
Harrison and Scott [2]. We shall further find, in analyzing the Higgs sector in the case of 
cyclic permutation symmetry, that the Higgs mass matrices are also circulants, making it 
easy to diagonalize them explicitly. 
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V. STRUCTURE OF THE HIGGS SECTOR 



We turn now to an analysis of the properties of the discrete chiral invariant Higgs po- 
tential of Eqs. (7a-c). We shall assume CP invariance and exact S3 cyclic permutation 
symmetry; when needed, we can take into account small deviations from these assumptions 
by adding perturbations to the locations of the Higgs minima. We begin our discussion with 
the three Higgs model, in which only the discrete chiral triplet (f) is present. Omitting the 
subscript (p on the coefficients, we have 

33 3 

-Cffiggs potential = (0n0n " V^ f - /^l X) ^n^n^l+l^n+l " X! l^n^n+lP 

n=l n=l n=l 

3 

-a X Re(j)'l(j)n+i(pi4>n-i ■ (11a) 

n=l 

Necessary conditions for this potential to be bounded below are evidently 

A>0, A-/Xi-/X2-q;>0 . (lib) 

Imposing the condition 

112 + a>0 (12a) 

insures that the Higgs potential is minimized when the three doublets all have the same 
form 

0n = f I , (12b) 

for a suitable choice of SU{2) gauge, with the consequence that one electroweak gluon (the 
photon) remains massless. Imposing the additional condition 

a>0 (12c) 

then forces the complex phases of the three expectations Q„ to be equal (up to discrete 
chiral rephasings) at the minimum of the potential; by a choice of U{1) gauge the overall 
common phase can be rotated to zero, and so the potential of Eq. (11a) is minimized at 
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(13a) 



with n given by 

n' = - — . (13b) 

A — fii — IJ,2 — a 

This minimum is not unique; because the potential of Eq. (11a) is invariant under the 
discrete chiral transformation of Eq. (3), equivalent minima are located at 



with c<Ji,2,3 any three distinct cube roots of unity, which can always be obtained by permuta- 
tion from the set a;, a;, 1. Despite the appearance of complex phases in Eq. (13c), there is no 
breakdown of CP invariance, because these phases can always be eliminated by the discrete 
chiral transformation that returns to the minimum of Eq. (13a). 

We note that although the potential of Eq. (11a) is similar in form to that studied by 
Bigi and Sanda [9], they choose o; < 0, in which case there are nontrivial relative phases 
(that are not just discrete chiral rephasings) between the three expectations 0.1^2,3 at the 
potential minimum, and CP is spontaneously broken. This case is not useful for our model 
building because numerical analysis shows that it leads to a mass matrix with one heavy 
family, and two other hghter families of equal mass. We shall make use of the possibility [10] 
of CP violation in multi-Higgs systems only in the context of the six doublet model, to be 
discussed shortly. 

To complete our discussion of the three doublet model, we must determine the Higgs 
masses. Expanding around the minimum of Eqs. (13a, b) to second order by substituting 



= 1,2,3 



(13c) 




(14a) 



into Eq. (11a), we find 



•Higgs potential 



Vo + V26 + V2e 



V2S 



V2e 




n. 



(14b) 
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A simple calculation shows that the matrices A,B,D are all circulants of the form 



A = Circ_(a^, b^, h^) , 

L> = Circ^(a^,6^,0 , 

with a^'^'^ and h^'^'^ given in terms of the Lagrangian parameters by 

= (2A - Ati - At2)2n2 - 2\v^ , = -(/xi + /xs + 2a)0^ 
^ 2{X - - 2Xv\ b^^-{n2 + ct)n^ , 

a^^{X-\a)n\ 6^ = -^(2/xi + 2/X2 + a)0^ . 



(15a) 



(15b) 



Since these matrices are all diagonalized by transformations based on the cube roots of unity, 
it is useful to introduce new bases defined as follows, 



02 



W 



w 



w 



0(2) 

V5(3); 

:(2) 



V6(3)y 



0(2) 

V0(^)y 
/5(i)\ 

V5(3); 

/6(1)\ 

.(2) 



-1 



V6(3)y 



'2 

V03y 

^2 
62 



(16a) 



with 



1 

/3 



/a; 


a; 


1\ 






a; 


1\ 


a; 


a; 


1 




a; 


CO 


1 


u 


1 




u 


1 


1/ 



iy^Circ_(a,6, c)!^ 



/ a + a;6 + a;c 


V 









a + cub + cue 

a + b + cj 
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WCirc^{a, b, c)W 







a + cub + cue 



a + cub + cue 

In terms of the new bases, Eq. (14a) becomes 



V 



\ 


a + b + c J 



(16b) 



(n) 



, n=l,2, 



(16c) 



and 



^^^^^ 



(16d) 



Substituting Eq. (16b) into both Eq. (14b) and the Higgs kinetic energy, and using 
Eq. (13b), we find for the terms quadratic in 5„ 



n=l,2,3 ^ 



= - E + («^ + 26'')^|5('f + (a^-6^)^(|5«|' + |5('f) 

n=l,2,3 2 2 2 

= - E + |(/^2 + «)Q^(|5Wr+|5(^f) . (17a) 

n=l,2,3 ^ ^ 

Prom Eq. (17a) we see that S^^^ is a charged massless Goldstone boson (which is absorbed by 
the Higgs mechanism into the longitudinal parts of the charged intermediate bosons), while 
(5(^'^) are two charged Higgs boson fields (each containing a positive and a negative charge 
state) , with mass squared 3(//2 + Q;)Q^. Similarly, we find for the terms quadratic in e„ 

- E lKen\' + V2, 

n=l,2,3 ^ 

= _ ^ l|a^e(")|2 + (a^ + 26^)-|e(3)|2 + (a^-6^)-(|e«|2 + |e(2)|2) 



n=l,2,3 



Defining new linear combinations e^^^ by 

eW = -^(eW±e(^)) , 



(17b) 



(17c) 
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and splitting e^^^ 6^=*=^ into real and imaginary parts, 6^^'=*=^ = ^i?'^ , Eq. (17b) takes 

the form 

n=l,2,3 ^ 

+(4A + 2/., + + la)n^l[{e^^^r + {e^f^n + + (4^^)^] ■ m) 

We see that e^^^ is a neutral massless Goldstone boson (which is absorbed by the Higgs 
mechanism into the longitudinal part of the neutral intermediate boson), while both 
and and both and e^^"*, are neutral Higgs states, with respective squared masses 
4Av^, (4A + 2//1 + 2/^2 + |q;)Q^, and |q;0^. Thus, the twelve states contained in the original 
triplet of Higgs doublets are accounted for as one neutral and two charged Goldstone modes, 
four charged Higgs bosons, and five neutral Higgs bosons. This information is summarized 
in Table I, which also gives the couplings of the Higgs bosons to fermions worked out in Sec. 
VI. 

We turn next to the properties of the Higgs sector of the six doublet model. Although we 
shall focus here on analytic results, we have also made numerical studies of the minima of 
the six (and three) doublet potentials, using the formulas and method given in Appendix A. 
Let us begin by assuming that the potentials Vi(0,?7) and ^2(0,77) of Eq. (7c), that couple 
the (j) and rj Higgs discrete chiral triplets, are very small. Then the minima of the Higgs 
potential are obtained by examining the degenerate minima of and V^, as analyzed in 
the three Higgs discussion above, and selecting those for which Vi + V2 is smallest. By a 
simultaneous Zq rephasing of and 77, we can always make the minimizing values of have 
the form of Eqs. (12b) and (13a), with Q and the coefficients A, v, a, //i, //2 in Eq. (13b) 
now carrying the subscript (j) to differentiate them from the similar formulas that hold for 
the Higgs field r}. There are now two distinct possibilities, depending on the values of the 
coefficients in Vi and V2. Suppose, for example, that all of the coefficients 7„ , C^mn in 
Eq. (7c) are negative; then Vi + V2 is clearly minimized if the expectations of rjn are all 
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relatively real to one another and to the expectations of 0„, that is, if 

?7n=^^j , (18a) 

with 

Ai = A2 = A3 = Q^ , (18b) 

with fljj given by Eq. (13b) with subscripts rj on all quantities. Suppose, however, that the 
coefficients in Vi and V2 are all positive; then the sum Vi + V2 will be made lower if we pick 
one of the degenerate minima of of the form of Eq. (13c), for example 

Ai = ujQjj , A2 = cuflri , A3 = . (18c) 

More generally, the necessary condition for Eq. (18c) to be a lower minimum than Eq. (18b), 
in the limit of small coupling of r] to 4>, is that Vi + V2 be smaller at Eq. (18c) than 
at Eq. (18b). Assuming exact cyclic permutation symmetry, which makes the following 
formulas independent of the value of the free index m, we find 

Vi'^^ ^ — 3f)^f2^ ^(C'lmn + C2mn + C'smn + C4rnn + C'smn + C'emn) ; 

n 

V2 ^ ^ ^ — 3^jn^4>^Ti + 3 y~^^[{Cjjnn + Csmri)^^^r] + {C^9mn + C'iOmn)^^^</) 

n 

+ (Cllmn + Ci2mn)^'l^'^] ' (19a) 

and 

^Eq.(18c) ^ 3^2^2^ , 

C = ^y^,\Cijnn ~ ^^C^mn + C'4mn)]^ + C'2mm + C'smm+l + C'emm+l 
^Eq.(18c)_Q . (19b) 

Thus, the necessary condition for Eq. (18c) to be the minimum is that 

^Eq.(18c) ^ ^Eq.(18b) _^ ^Eq.(18b) _ ^^g^-j 



We shall henceforth assume that Eq. (19c) is satisfied; as already noted, this is automatic 
in the case when all of the coefficients in Vi and V2 are positive, but the general condition is 
much less restrictive, requiring only that the coefficients lie on one side of a hyperplane in the 
space of Vi,2 coefficients. When Eq. (19c) is satisfied, CP invariance is spontaneously broken 
through the rj Higgs expectations, and we shall see in the next section that simultaneously, 
the r] expectations have the correct form to generate nonzero second family masses. 

Let us next consider what happens when Vi and V2 are not infinitesimally small. Still 
maintaining cyclic permutation invariance, let us first consider the case in which Vi is large, 
but V2 remains nearly zero. Then from the formulas of Appendix A, we find that the 
derivatives of the potential vanish when one assumes Eq. (13a) for the expectations (with 
of course replaced by Q^) and Eq. (18c) for the rj expectations, for suitable minimizing 
values of and fl^j. In other words, we find the correct minimum by first substituting 
Eqs. (13a) and (18c) into the Higgs potential, and then minimizing the resulting simplified 
expression with respect to and Q^. Substituting Eqs. (13a) and (18c) into Eq. (7a) gives 

^^Cmggs potential = A^Q^ - 25^1^2 ^ ^^^4 _ 2^^q2 ^ c'fijfij + coustaut , (20a) 

with C given in Eq. (19b), and with the remaining coefficients given by 

A^^X^- - ii2^ - a^, = A^v|, ^^(t),v ■ (20b) 

Minimizing Eq. (20a) with respect to Q^, fl^ gives a pair of simultaneous linear equations, 
with the solution 

^ ~ A^Ajj — |( 

In order for both and rj to develop nonzero vacuum expectation values we must have 
> 0, > 0, which in the case when the denominator in Eq. (20c) is positive requires 
that C be restricted by 

-2AIaI <C <2Mm{^B^,^B^,AlAl) . (20d) 
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Because Vi is invariant under independent overall phase rotations of and rj, in the limit 
when V2 is strictly zero the minimum of Eqs. (13a) and (18c) is part of a one parameter 
U{1) family of equivalent minima, of the form 

(Oi,02,^^3) = (l,l,l)f^<A , 

(Ai, A2, A3) = (oJ, ou, 1)Q^ exp{i0) , (21) 

with the angle 9 arbitrary. When V2 is nonzero but very small, the U{1) degeneracy with 
respect to 6 is broken, and the minimum has the form of Eq. (21) with a definite value of 
9 determined by the Higgs Lagrangian parameters. A perturbative analysis in powers of 
V2 shows that to first order in V2 the degeneracy in 9 is unbroken [because the final line of 
Eq. (19b) remains valid for general 9], but that at second order in V2 a nontrivial condition on 
9 is obtained and the degeneracy is broken. Numerical minimization of the Higgs potential, 
using the method of Appendix A, shows that general values of 9 can be attained at the 
minimum for generic Lagrangian parameters. As V2 increases, there are relative phase and 
small magnitude corrections to the minima of Eq. (21); when the assumption of cyclic 
permutation symmetry is relaxed, these magnitude corrections become more pronounced. 

To conclude our discussion of the six Higgs model, let us discuss the Higgs mass spectrum, 
assuming both exact cyclic permutation symmetry and the weak coupling limit in which both 
Vi and V2 are very small. We parameterize the expansion of 0„ and rjn around the minimum 
as 



(22a) 



Vn^\ ^ ] exp{i9){uj,uj,l)n , (22b) 

where we have used the notation {x, y, z)n to indicate x for n = 1, y for n — 2, and z for 
n — ?>. Because the overall phase 9 and the discrete chiral phases (oJ, a;, 1)„ drop out of VJ,, 
for the non-Goldstone modes we get simply two copies of the nonzero mass modes found 
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in Eqs. (17a-d) in the three Higgs case, apart from adding subscripts or superscripts 0, 77 
to distinguish the (f> and rj sectors, as summarized in Table II. In computing the Yukawa 
couphngs of the rjn Higgs modes, the phases in Eq. (22b) play a role. Making transformations 
analogous to Eqs. (16a) in the three Higgs case, with ^ in the following formulas either or 
77, we have 



f4\ 
\4j 



= w 



w 



w 



.(2) 



^(2) 



w- 



w- 



/6\ 

6 

4 
4 

\4/ 



In terms of the new bases, Eq. (22a) becomes 
and 



J_An) 



n 



1,2, 



while taking into account the extra phases, Eq. (22b) becomes 



(3) 



V2 ^ 



exp(i^) 



(23) 



(24a) 



(24b) 



(24c) 



■>(!) 



and 



^^(2) 
V2 »? 



^^(3) 



exp(i6') 



exp(i^) 



(24d) 



(24e) 
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The fact that appears in 77^^) rather than in 77^^) is directly related, as we shall see in the 
next section, to the role of the f] Higgs bosons in giving rise to second family masses. 

For the Goldstone modes, the situation is more complicated, because the 4> and r] sectors 
interact even in the weak coupling limit. If V2 were exactly zero, as noted above we would 
have an extra U (1) symmetry, and we would get two copies of the Goldstone modes as well. 
But for nonzero V2 this U{1) degeneracy is broken, and we are left with just one set of 
Goldstone modes, corresponding to the remaining invariance of the Higgs potential under 
simultaneous overall rephasing of 0, rj, while the three Goldstone modes related to the relative 
phase 6* of and r] become massive pseudo-Goldstone modes, with squared masses that are 
proportional to the magnitude of V2. The decomposition of 5'^^ and ef^^^ into Goldstone 
and pseudo-Goldstone modes is made unique by the facts that (i) these represent orthogonal 
degrees of freedom, that are simply rotations from the original modes 5^ ;^ and e},/,^, and 
(ii) the Goldstone modes correspond precisely to a uniform infinitesimal phase rotation of 
(f), rj, which specifies the infinitesimal modes to which the pseudo-Goldstone modes must be 
orthogonalized. Since the expectations of 4>,r] may have unequal magnitudes ^2,^,^2^, we 
see from Eqs. (24a-d) that an overall infinitesimal phase rotation makes a contribution to 
5^^^ that is flri/^(j) times as large as the corresponding contribution to 5'^\ and similarly 
makes a contribution to e}^ that is VLr^/Vt^ times as large as the corresponding contribution 
to ef^. We thus find, denoting the Goldstone and pseudo-Goldstone modes respectively by 
the subscripts G and PG, and as before using the subscript I to denote the imaginary part. 



O -U O ^(3) o ^(3) Q (3) 

The corresponding quadratic terms in the Lagrangian are 

(25b) 
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The perturbative contribution to the pseudo-Goldstone masses, relative to the Higgs masses 
calculated above, will have the general magnitude (suppressing all subscripts) 

We have not attempted to calculate explicit perturbative formulas for the pseudo-Goldstone 
masses, both because these will be rather complicated given the complexity of V2 and be- 
cause, as argued by Weinberg [11], there are likely to be significant nonperturbative correc- 
tions of order gMw, with g the electroweak gauge coupling and M]y the electroweak boson 
mass. 



VI. HIGGS COUPLINGS AND MASS AND CKM MATRICES WHEN CYCLIC 

PERMUTATION SYMMETRY IS EXACT 

We proceed now to study the Yukawa couphngs of the Higgs fields, and the mass ma- 
trices generated by their vacuum expectation values, when cyclic permutation symmetry is 
exact. Thus, in this section we shall assume that the Higgs potentials have the cyclically 
symmetric form analyzed in detail in Sec. V, and we shall take the asymmetry parameters 
(^(im of ^'i^- (6b, c) to vanish. As a consequence, the 3x3 matrices of Eq- (6b) are all 
retrocirculants, and are independent of the labels ^, /, 

P/i-Circ^(0,l,0) , 
P/2 = Circ^ (0,0,1) , 

P/3 = Circ^(l,0,0) . (26) 

Substituting Eq. (23) for ^1,2,3, with ^ = 4>,r], into the first line of Eq. (6a), we get for 
f ^d,e, 

- 9i{Pl^''<t>^'^ + Pl^'U^'^ + Pl^'U^'^) + giiPi^H'^ + P^^'V'^ + P^^V^) . (27a) 
Here we have defined 
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= w 


pf 




pf 


Wv 











(27b) 



with W and W ^ as given in Eq. (16b). Defining CP conjugates of by 
the CP conjugate of the first group of equations in Eq. (23) is 



(28a) 

















= w-^ 




|(2) 


= w 




Via J 










U~3/ 



(28b) 



''1,2,3, Vi,2,3 in the second Une of Eq. (6a), we get for f = u,^, 



Using this for (pi 

$/ = 4(i^(^)<^(i) + Pl^'^cp^^ + pI^4>^^) + gliP^^'^'^ + Pi^'¥'^ + Pi^'¥'^) . (28c) 

Substituting the retrocirculant forms of Eq. (26) into Eq. (27b), we can write the matrices 
P/*^^'^'^^ as retrocirculants. 



,/{2) 



j/(3) 



f^') = -^Circ^(l,a;,a;) , 



V3 
1 



Circ^(l, a;, a;) , 



^ -^Circ^(l,l,l) . 



(29a) 



Let us now use Eq. (10b), which asserts that p^^^'^'^^ are all diagonahzed by the same 
bi-unitary transformation constructed using Ul, Uj^ of Eq. (10a), 



ULP^Ui = V3diag(l, 0, 0) = V3M(i) ^ 
UlPI^'^^uI = v/3diag(0, 1,0) = ^SM^^) , 



UlPFuI = V3diag(0,0, 1) = ^3M^^^ 



(29b) 



Clearly, the natural thing to do now is to rotate to new fermion bases using the same matrices 
Ul, Ur, by introducing primed bases defined by 
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fn = Uy'ji, f = d,u,e,iy . (30a) 

Since the fermion kinetic energy of Eq. (4c) does not couple left to right chiral components, it 
has the same form in terms of the primed bases as in terms of the original ones. Substituting 
Eqs. (27a), (28c), and (29b) into the Yukawa Lagrangian of Eq. (5b), we get finally 

>CYukawa = Ql^'^^A + Ql*"^^^ + + V^L*'^^ + ^djoint , (30b) 

with the 3x3 matrices defined by 

*^ = Ex/3(4'^(^)+5,VW^ f-d,e , 

*^ = EV3(40^'^+4r^('^)MW, f^u,u . (30c) 

e=i 

On substituting Eqs. (24a-d) into Eq. (30c), we can read off both the mass matrices 
and the Yukawa couplings of the physical Higgs states. The mass matrices are obtained by 
keeping only the vacuum expectations of r]^^\ that is, by setting 

0(1.^) ^ 0, ^ 0, - f ° V 0^^)- f"^^^") , (31a) 

and 

/ \ / ^/3n„exp{-ie)\ 

^(1,3) ^ 0, fj^'''^ ^ 0, 77(2) ^ , 77^^) ^ , (31b) 

Vy3Q^exp(ie)/ V / 

giving 

>Cn.a.s = 5;(3^^Q,exp(i^)M(2) + 3gp^M^^^)d'j, 

+«i(3^;;i],exp(-z^)M(2) + 3g;n^M^^^)u'j, 
+ei(3^^Q,exp(i^)M(2) + 3g;n^M^^^)e'j, 

+i7]^(3^^f]^exp(-i^)M(2) +3^^f)<^M(2))i/R + adjoint . (32a) 

Identifying M*^^'^'^^ respectively as the projectors on the first, second, and third family states 
in the primed basis, we read off from Eq. (32a) the masses 
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Mt = ?,gl%, Me = 3^^Q„ M„ = , 
M, = 3gp^, Ms = Sgprj, Ma = , 
Mr = Sg;n^, M^ = Sg'fy, Me = , 

M,^ = 3g;n^, M,^ = 3^;;fi„ M,^ = . (32b) 

We see that in the three Higgs model, only the third family gets masses, with the first 
two families remaining massless. The same is true in the CP conserving phase of the six 
Higgs model, in which the 77 expectations are given by Eq. (18b) rather than Eq. (18c); 
in this phase, the projectors M^^^ in Eq. (32a) are replaced by projectors M^'^\ and the 77 
expectations simply make additional contributions to the third family masses. On the other 
hand, in the phase of the six Higgs model that spontaneously violates CP as in Eqs. (18c) 
and (21), the factors oJ, a;, 1 in Eq. (18c) give rise to the projector M^^^ for the second family 
states, which then receive masses. The hierarchy between the masses of the second and third 
family charged leptons is attributed, in the six Higgs model, to a systematic tendency of the 
77 Higgs bosons to have smaller Yukawa couphngs to the charged fermions than those of the 
(j) Higgs bosons. 

To get a feeling for the magnitudes involved, we note that the Higgs boson expectations 
generate mass terms for the gauge bosons given by 

with 

^^ = 2^(|(0„)|2 + |(r7.)r)=6(nJ + Oj) . (33b) 

n=l 

Empirically, v ~ 247GeV; assuming, as we shall in the fits below, that Q,^ and fi^, are 
approximately equal, we then find ~ ~ 71GeV. The Yukawa couplings needed to 
reproduce the observed charged fermion masses are then given in the six Higgs model by 

gl ~ 0.81, ~ 0.0061 , 
gl ~ 0.020, ~ 0.00094 , 

gl ~ 0.0083, gl ~ 0.00050 . (34a) 
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In the three Higgs model, is a factor y^2 larger than in the six Higgs model, and the 
Yukawa couplings are correspondingly a factor -^2 smaller than in Eq. (34a), 



As we have seen, because the mass matrices in the cyclically symmetric limit are retro- 
circulants, we were able to diagonalize them with universal, flavor independent matrices 
Ul, Ur. This has the important consequence that when cyclic symmetry is assumed as a 
leading approximation, the corresponding approximation to the CKM mixing matrix is unity, 
a welcome feature since the observed CKM matrix is close to unity. A related welcome fea- 
ture of the cyclic approximation is that there are no flavor changing neutral currents, which 
again accords with the fact that these are observed to be highly suppressed. To obtain 
realistic non-unit values for the CKM matrix, we will have to go beyond the cyclic approx- 
imation by including nonzero asymmetries /^^^^ as in Eq. (6b), but we shall then also have 
to estimate the magnitude of the flavor changing neutral current effects produced by these 
asymmetries. This will be the agenda of the next three sections. 

Before proceeding with this analysis, we note that the leading cyclic approximation to the 
Yukawa couphngs of the physical Higgs bosons can be read off from Eqs. (30b, c) together 
with Eqs. (24a-d) and (25a). We see that 



gl ~ 0.57, 4 ~ 0.014, ~ 0.0059 



(34b) 



a(1) A(2) .(2) 

a(2) .(2) X(3) (3) 
.(3) r(l) (1) 



couple only to the first family 



couple only to the third family 



couple only to the second family 



(35a) 



which by Eq. (17c) imply that 




couple only to the first and second families 



eljji j couple only to the first and third families 



(35b) 



and by Eq. (25a) imply that 



5pQ, CpQ couple only to the second and third families 



(35c) 
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The presence of 21 Higgs bosons in the six Higgs doublet model [eight charged Higgs bosons 
S^J^\ ten neutral Higgs bosons e^],^, e^^^, and C/^i*^, plus two charged pseudo-Goldstone 
Higgs bosons 6pQ, and one neutral pseudo-Goldstone Higgs boson Spq], together with the 
pattern of predominant fermionic couplings given in Eqs. (35a-c) and summarized in Table 
n, are a distinguishing feature of the model that should be testable in experiments at the 
next generation of accelerators. 

VII. FIRST ORDER BREAKING OF CYCLIC PERMUTATION SYMMETRY 

We now set up a perturbative scheme to study the effects of the breaking of cyclic 
permutation symmetry. In the three Higgs model, we will also allow CP noninvariance of 
the Lagrangian, by allowing the phases ip in Eq. (7b) to be nonzero, and by allowing the 
Yukawa couplings to be complex. In the six Higgs model, we will impose CP invariance on 
the Lagrangian, but will work in the phase that spontaneously breaks CP. Two types of first 
order small corrections will be introduced. The first are corrections to the Higgs vacuum 
expectations, arising from a lack of cyclic symmetry in the Higgs potential. In the three 
Higgs model, this results in replacing Eq. (13a) by 

fi„ = l](l + 5„), n = l,2,3 , (36a) 

where the Sn are small corrections that can be complex, and where we impose the condition 

E^n = (36b) 

n 

to avoid duplicating information contained in the overall factor and the overall phase that 

has been eliminated by a gauge transformation. In the six Higgs model, we have analogous 
corrections to the first line in Eq. (21), 

n^^n^il + Sn), n = 1,2,3, J2^n^o , (37) 

n 

where the 5„ can again be complex when the potentials Vi,V2 that couple (p to rj are not 
neglected. In principle, there are also asymmetry corrections to the second line of Eq. (21), 
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which gives the r) expectations. But these are always suppressed by a factor gf^/g^, which 
according to Eq. (34a) is at most of order 0.06, and so will be neglected in what follows; 
that is, we treat grj/gtf, here as if it were also a first order small quantity. The second type of 
first order small corrections are the asymmetry parameters P^^^ ^^S- (6b, c), which are 
complex in the three Higgs model when explicit CP violation is permitted, but are real in 
the six Higgs model when CP invariance is imposed on the Lagrangian. Again, in principle 
there are analogous asymmetry parameters for the rj Yukawa couplings, but the effect 
of these is again suppressed by a factor g:lj/g^ and so they will be neglected. This itemization 
of corrections defines the model that we shall study in first order perturbation theory. 

Since the zeroth order problem, that was analyzed in Sec. VI, is brought to diagonal form 
by the bi-unitary transformations of Eqs. (29b) and (30a) based on the matrices Ul,Ur of 
Eq. (10a), we shall make this transformation at the outset. In the primed fermion basis, 
the zeroth order mass matrix is still given by Eq. (32a), but now there will be first order 
corrections from the 5's and /5's introduced above. Since we are regarding g^/gjj, as effectively 
a first order correction, it is convenient to group it with the other first order terms. Starting 
again from Eqs. (23), (27a, b), and (28c), we then find for the extension of Eq. (32a) to 
include all first order corrections, 

Cr.^^ E 7l4^^<^(3M(^)+(j^)/;, , (38a) 

f=d,u,e,i/ 

with cr-'^ a 3 X 3 matrix with matrix elements given by 

(^11 = ^^11 +4 + ^^2 + ^(^1 , 

^2 = + + + ^4 + ^4 , 

4 = , 

<^/m = ^/^L> 7^ ^ ■ (38b) 

The further quantities appearing in Eq. (38b) are defined as follows. The quantities are 
given, in terms of the Sn introduced in Eqs. (36) and (37), by 



29 



S( = dl, 6^ = 6;, di = d;, f = u,u . (39a) 
The quantities are defined by 

witli tlie + sign holding for f = d,e and the — sign holding for f = u,^. Finally, the A*{^'s, 
when multiplied by the factor of 1/3 in Eq. (38b), are the asymmetries /?^^^ reexpressed in 
the primed fermion basis; suppressing the subscript on the /3's, they are given by 

/^fi = pL + pL + Pi2 + ^{^(2 + pL + pL) + coipL + pL + pL) , 
d2 = P(i + Pis + Pi2 + ^{Pi2 + pL + pL) + ^{pL + pi2 + pL) , 



1^(2 = p(i + pL + pL + ^{Pi2 + pL + pL) + ^{pL + Pi2 + 

/^^l = Pil + /324 + PL + ^(/5f2 + + PL) + + Pi2 + /^fa) 

/^fs = Pil + /3f2 + /S/a + HPL + + PL) + + /9|2 + PL) 

ds = + P(2 + pL + ^{pL + pL + PL) + ^{pL + pL + pL) 



f^ii = Pil + PL + + ^(/9f2 + PL + PL) + u^iPis + /Sis + /Sla) , 

^2 = + pL + pL + ^{P(2 + PL + pL) + ^{pL + pL + pL) ■ (39c) 

We remark that since CP invariance requires the /3's to be real, the condition for CP invari- 

ance, when expressed directly in terms of the //'s, is //{* = //22, A*i2 — A*2i) Mia — A*23) 

f* f 

1^31 = A''32- 

Defining 

M'f = 3M(3) + , (40a) 
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we must now find tfie bi-unitary transformation matrices U[,Uj^ for which U[M'j^Uj^ is 
diagonal, with the eigenvalues ordered in absolute value, for each flavor /. The fermion 
basis states that are mass eigenstates are then related to the primed basis by 



•mass 



fR-U{^fr^ f^d,u,e,u , (40b) 
and the CKM matrix Uqkm is given as usual by 

Uc^M^U^Ut . (40c) 

We shall now develop a perturbative procedure for calculating C// The first observation 
to be made is that wc arc deahng with a degenerate perturbation problem, since the zeroth 
order mass matrix 3M = 3diag(0, 0, 1) has eigenvalues for the first two primed basis 
states. As a consequence, the 2x2 submatrix of C// spanned by these states is zeroth 
order in the perturbation o"-^, with only the off-diagonal elements coupling to the third basis 
state of first order. Thus, we find a natural reason in our model why the CKM mixings of 
the first and second family states should be larger than the mixings of the first and second 
families with the third family. 

We shall deal with the zeroth order 2x2 submatrix by calculating it exactly. Let Vl j^ 
be the 2x2 matrices that bring the 2x2 submatrix of to diagonal form, 

Vl{''' '] . (41a) 

with the magnitudes of the eigenvalues ordered as | k{ | < | 1 • The explicit construction of 
^L,R is given in Appendix B. It is then straightforward to show that to first order in small 
quantities, C//^^ are given by 



a 



( V[ -W[\ 



L 3 "L 

„f \ t 
^13 



(41b) 
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and 



Defining 
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/t 

R 



/* \ t 



/* 



/t 



R 



'32 



UlM'^U^J 



f k{ 0\ 


V 3v 



(41c) 



(41d) 



(42a) 



the corresponding first order accurate expression for the CKM matrix is given by 



U, 



CKM 



VcKM 



a 



KM 
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d \ t 
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'23 , 



a 



13 



'23 , 



a. 



13 



(42b) 



23 , 



'23 . 



VcKM 1 

Although Eq. (42b) is useful for analytic study of the CKM matrix, in our numerical work 
we shall simply compute directly from the definition of Eq. (40c). We shall also, in the 
numerical work, use slightly more accurate forms for t// in which only the square of the 
ratio of second to third family masses (|ac2|/3)^ is assumed to be small; the relevant formulas 
are given in Appendix C. 



VIII. HIGGS EXCHANGE CONTRIBUTIONS TO THE Kl - Kg MASS 

DIFFERENCE 

As pointed out in Sec. VI, when cyclic symmetry is exact, Higgs boson exchange in 
our models does not produce strangeness changing neutral current effects. However, once 
we include cyclic asymmetries, such effects become possible and we must be sure that their 
magnitude does not exceed known experimental limits. Since, in the context of extensions of 
the Higgs sector, the most stringent bound on strangeness changing neutral current processes 
comes [13] from the second order weak — Kg mass difference, we shall consider only 
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this process, and shall calculate the contribution to its matrix element arising from Higgs 
exchange within the perturbative framework set up in Sec. VII. 

We saw there that, because the zeroth order mass matrix is degenerate in the subspace 
spanned by the first two families, the mixing matrices within this subspace are zeroth order 
rather than first order in the perturbation, and therefore strangeness changing neutral cur- 
rent effects can already appear at zeroth order in perturbation theory. What we shall do in 
this section is to calculate this zeroth order contribution to the Kl — Ks mass difference, 
neglecting all terms of first and higher order in the asymmetric perturbation. Our starting 
point is thus the Yukawa Lagrangian of Eqs. (30b, c) in the primed basis, of which the term 
relevant to mixing of the d and s quarks is 

QlEV3(4'/'^'^ + *^'^)^^'^4 + adjoint . (43a) 

Substituting Eq. (40b) relating the primed to the mass eigenstate bases and using the ap- 
proximation of Eqs. (41b, c) for Vl j^] also substituting Eqs. (24a, b) for the (p^^^ and keeping 
only the neutral Higgs pieces; and finally also neglecting terms of first and higher order in 
the asymmetric perturbation, we get the effective Lagrangian 

>Cscnc ^ 3r^^4E 4'^^lMS2^«'^^"^ + adjoin . (43b) 

In Eq. (43b), the subscript 2 x 2 on the projectors indicates their restriction to the subspace 
spanned by the first two families, and the column vector d will be understood to have been 
truncated from three to two components, corresponding to the first two families. Finally, 
reexpressing e^'^^ in terms of the modes e^"*^^ defined in Eq. (17c), sphtting these into real 
and imaginary parts, and explicitly including the adjoint term (our 7 matrix conventions 
are 75 = 7^, 7° = 7°+, (7°)^ = 1), we get 

+5"^"^^<{[(6« + eg) - .(eW + e[',^)nM^^iX' 
+[(4""^ - 4«) - ^(4/^ - ^^^WnM^i^Vt'Ml - l,)d--' . (43c) 
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To facilitate the remaining calculation, it is convenient to rewrite Eq. (43c) in the form 

(A^) + 5i+)75)(i"^-^ + d^'^e^,) (A^ ) + 5}-)75)ci'"-^ . (44a) 

Using the facts that 

MS2 + M^i, = 1^ ^ ° j = 1, MiJ^ - MiJ^ = ^ ^ °i ) = ' ^^^^^ 
we find that the 2x2 matrices A^^J, B^^j appearing in Eq. (44a) are given by 

(44c) 

A'' = ^9lWtp,V'^ - ^gt^Vip,Vt\ 4"^ = ^9Kp.y'^ + ^9t^yiP^yt' ■ (44d) 

Letting d and s denote, respectively, the down and strange quark eigenstates, the two com- 
ponent column vector (p^^'^ has the structure 

^mass , (45a) 

and so for any 2x2 matrix iV, we have 

d'^^'NcT^ = dNnd + dNi2S + sN2id + SN22S . (45b) 

Hence the strangeness changing terms of Eq. (44a) involve only the 12 and 21 matrix elements 
of the matrices in Eq. (44c), and can be compactly written as 

^tlf = E E [detl{A% + Bf,2l.)s + sefAA% + 4l75)ci] . (45c) 

p=± F=R,I 

So for the amphtude T for the AS" = 2 process s + s— >d + dwe find, summing over the 
exchanges of Higgs eigenmodes with squared masses M'§-^^ , the formula (valid up to an 
overall phase) 
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r = E E ^(4^2 + 4^1^5)^-4)^(412 + ^iWs)^ , (46a) 

p=±F=R,I Mp 

while from Sec. V and Table II we find for the squared masses, 
M^W ^ ^2(-) ^ ^^^^ ^ 2;,!^ + 2//2^ + ^q;^)Q| , 

mT' - Mf-^ ^la,^l nl^- ^-fi . (46b) 

The Higgs exchange matrix element {K\T\K)iiiggs for the K ^ K transition, in the vacuum 
saturation approximation, then is given by 

|(K|r|F)Higgs|^A^|(i^|^75s|0)nL'Higgs| , (47a) 

with i^Higgs given by 

^Higgs = E E ^ , (4rb) 

p=±F=R,7 -'"F 

and with iV = 8/3 a Wick contraction and color factor. 

We wish now to compare the amplitude of Eq. (47a) with the intermediate boson loop 
diagram contribution to the Kl — Kg mass difference calculated by Gaillard and Lee [14], 
which in vacuum saturation approximation is in satisfactory agreement with experiment. 
The Gaillard and Lee result is 

|(i^|r|Z)GL|=^iV|(i^|d7^75s|0)nL>GL| , (48a) 

with \Dqi\ given by 

I^glI = 0M,2s?2 , (48b) 

in terms of the Fermi constant Gp-, the charm quark mass Mg, and the sine of the Cabibbo 
angle Si2 = smOc- To compare Eq. (47a) to Eq. (48a), we need the ratio of the pseudoscalar 
current to the axial vector current kaon to vacuum matrix elements, which can be estimated 
by standard current algebra methods (see, e.g., Shuryak [15]) to be 

'W75.|0)P ^ (OMO)! ^ ^M^^2 ^ 11 (49) 
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with Mk and fx the kaon mass and decay constant and with Mg the strange quark mass. 
Combining everything, we find that the condition for the Higgs exchange contribution to 
the Kl — Ks mass difference not to exceed the Gaillard and Lee estimate is 

A^mi ^ GeV^ ' 

which will be used as the strangeness changing neutral current constraint in the fits of the 
next section. 



IX. NUMERICAL FITS OF THE THREE AND SIX HIGGS MODELS TO THE 

EXPERIMENTAL DATA 

In order to fit the models to the experimental data, we follow the standard procedure of 
minimizing a "cost function" C, constructed as follows, 

C — Cmass + C'cKM + C'scnc + ^'parameter ) (51a) 

with the pieces referring respectively to the constraints placed by fitting the masses, fitting 
the CKM mixing angles, obeying the strangeness changing neutral current bound of Eq. (50), 
and keeping the asymmetry parameters as small as possible. Before giving further details, 
we describe the general search method employed. We perform all fits using the minimization 
routine "powell" of Press et. al. [12] . As given in [12] , this routine works well for the six Higgs 
model where the degeneracy between the first and second families is already broken, before 
inclusion of the asymmetry parameters, by the r] Higgs couplings. However, in the three 
Higgs model, it is not a priori specified which states become the first and which become 
the second families, and so eigenvalue crossings can occur in the course of the iteration 
which result in discontinuous behavior of the cost function. This causes a problem with 
the bracketing routine "mnbrak" of [12], which occasionally gets stuck in an indefinite loop. 
The fix is simply putting an iteration counter into "mnbrak" , to force an exit with a default 
bracketing (specifically, in terms of the quantities defined in "mnbrak" , c — a, fc — fa) if 
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convergence to a bracketing is not attained in A^max passes through the loop. We found the 
same results in the three Higgs model with A^max = 5 as with A^max — 30, indicating that a 
bracketing is attained very rapidly, or not at all. As an additional check, we verified that 
the original and the modified versions of "mnbrak" give identical results for the six Higgs 
model, where level crossings and associated discontinuous behavior do not occur. 

Let us now turn to the construction of the various cost function terms in Eq. (51a), 
working throughout in units where 1 GeV=unity. For the mass cost function, we use a 
standard chi squared function constructed from expected values of the masses and their 
estimated errors, including the electroweak mass parameter v of Eq. (33b). To prevent 
the chi squared for certain very accurately known masses (such as the electron mass) from 
dominating the fits, we truncate these masses to a few significant figures and use enlarged 
error estimates. In the six Higgs model we also add a term that favors fits with fi^ ~ il^, 
since this degeneracy plays a role in the extension to neutrino mixings discussed in the next 
section; in practice, we find that this term has very little effect on the fits, since nearly equal 
values of f2<^ and are favored even in its absence. (This term is omitted in the three 
Higgs model, where it is not relevant.) Adding these contributions, we have for the mass 

cost function Cmass; 

/M„- 0.005\2 /Mc-1.3\^ /Mt-173 



/ Mrf-0.01 \^ / M, -0.2 ^ / Mb -4.3 - ' 
V 0.005 ) \ 0.06 ) V 0.2 
/Me- 0.00051\2 /M^- 0.1057\2 /M^ - 1.777 



V 0.0001 ; V 0.001 ; v o.ooi 

2 



/ [6(Q^ + Q^)]^- 247. 



+ 

\2 



+(l^^-fi^)^ . (51b) 

For the strangeness changing neutral current cost function, we use a chi squared function 
with expectation zero and standard deviation equal to the bound of Eq. (50), 

C-c = (2.1 x^'lO-i^) ■ 
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To set up the CKM cost function, we make the standard rephasings to put the CKM 
matrix in the form 

/ 1 Si2 si3e-'^''\ 



CKM 



(52a) 



"•§12 1 S23 

V-Si3e^^" -S23 1 J 
to first order accuracy in small quantities, and then construct a chi squared function from the 
expected values and estimated errors for S12, S13, and S23- Although the CP violating angle 
5i3 has not been reliably determined experimentally, it appears likely that it is appreciable, 
so we also include a chi squared term requiring | sin^ial to be equal to 0.6 ± 0.3, giving 

'si2-0.221\2 /si3 - 0.0035\2 



y / 513 - 0.0035 Y 
J V 0.0009 / 
/S23 - 0.041\2 /|sin(5i3| -0.6~ 



2 



Altogether, then, there are 11 quantities to be fitted in Cmass, 1 to be fitted in Cgcnc, and 4 
to be fitted in Cckm, for a total of 16. 

Let us now count the numbers of parameters in the two models, and establish the cost 
functions for the parameters. Despite its increased complexity in terms of particle content, 
the six Higgs model has the smaller number of parameters, since it violates CP only spon- 
taneously and so all Yukawa couplings appearing in the Lagrangian are real. Altogether, 
there are 37 parameters that enter into the iterative fit for the six Higgs model. These are 
the 4> and 77 expectations il^ and flj^, the real parts of the Yukawa couplings g^^^, f = u,d, e, 
the complex asymmetry parameters 5i^2 introduced in Eq. (37), the angle 9 of Eqs. (21) 
and (39b), and the real asymmetry parameters /S^^n, f — u,d,e, m + n < 6 introduced in 
Eqs. (6b, c). The parameters A</,, /xi0 + /X20, and a^, which enter the calculation only through 
their appearance in the Higgs masses in the strangeness changing neutral current constraint 
[see Eqs. (46-47)], were fixed at the respective values 1, 0.3, and 0.3, and were not iterated. 
To construct the cost function for the iterated parameters, we note that no additional con- 
straint is needed for the expectations Q^,^ or the Yukawa couplings g^ ,^ because these are 
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already adequately controlled by Cmass of Eq. (51b). For the remaining parameters we use 
the cost function 

3^ 

r^Smn 



pajrameter 



E 

n = 1,2 



'^parameter 



+ E 

m, n m + n < 6 



'-'^parameter 



+ 



e 



6.28 



(53a) 



f = u,d,e 

In Eq. (53a) the exponent e and the width (Jparameter are parameters of the fitting procedure, 
which effectively set up a model for how the small asymmetries are distributed. We were 
able to get satisfactory fits for both e = 1 and e = 2, but convergence was much slower for 
the latter, suggesting that e = 1 is a model in closer correspondence to the experimental 
data, and we shall only present the e — 1 results in the discussion below. To initialize the 
six Higgs minimization search, we started from Vt^i = = 70.7, the values of Eq. (34a) for 
the Yukawa couplings g^^^^, zero for the asymmetry parameters and zero for 9. 

Because the three Higgs model, to give a CP violating CKM matrix, must violate CP 
explicitly, its Yukawa couplings and Yukawa asymmetries can have imaginary parts, and 
so there are 57 parameters that enter into the iterative fit. These are the expectation 
Jl<^, the real parts of the Yukawa couplings g^^, / = u,d,e and the imaginary part of 
(since gf^'^, which are not involved in the strangeness changing neutral current constraint, 
enter only through their absolute values, they can be rephased to be real), the complex 
asymmetry parameters 5i^2 introduced in Eq. (36a), and the complex asymmetry parameters 
f^imn, f = u,d,e,m + n<6 introduced in Eqs. (6b, c). Again, the parameters A^, + 
and a^i), which enter the calculation only through the strangeness changing neutral current 
constraint, were fixed at the respective values 1, 0.3, and 0.3. To construct the cost function 
for the iterated parameters, we note that again no additional constraint is needed for the 
expectation Q,^ or the real parts of the Yukawa couplings g^, because these are adequately 
controlled by Cmass of Eq. (51b). For the remaining parameters we use the cost function 
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pajrameter 



E 

n = 1,2 



d, 



'^parameter 



E 

m,n m + n < 6 



4>mnF 



'^parameter 



+ 



0.028 



(53b) 



F^R,I 

The width 0.028 governing g^j is chosen here as twice the natural magnitude of according 
to the estimate of Eq. (34b), so as to bound g^j but not overly restrict it, much as the 
width for 9 in the six Higgs model is chosen in Eq. (53a) as twice the maximum magnitude 
TT of 1^1 . Again, the exponent e and the width CTparameter are parameters that model how the 
small asymmetries are distributed. For comparison with the six Higgs model fits, we shall 
again only present e — 1 results in the discussion that follows. To initialize the three Higgs 
minimization search, we started from fl^ — 100, the values of Eq. (34b) for the real parts 
of the Yukawa couplings g^, zero for the imaginary part of g^, and zero for the complex 
asymmetry parameters (5„, 

We begin by presenting results for the six Higgs model. In any fitting procedure involving 
more parameters than quantities to be fit, one has to worry about overfitting, and we deal 
with this in the following way. As we shall see shortly, the most sensitive aspect of the 
fitting procedure for the six Higgs model is getting the CKM parameters correct, and so we 
take the cost function subcomponent Cckm as a measure of overfitting. Making a series of 
fits using the cost function of Eq. (53a) with e = 1, as a function of the width (Jparameter, 
we find that the value of Cckm is a monotone decreasing function of the width. For very 
small values of the width (i.e., asymmetries restricted to have very small values) we find a 
value of Cckm much larger than 4, the number of fitted CKM matrix degrees of freedom; 
for large values of the width we find values of Cckm much less than 4, indicating overfitting. 
We take as "good" fits ones resulting from widths ^parameter that yield a Cckm of order 4; an 
example of such a fit, with (Jparameter — 0.03, is given in Table HI. This fit, which was attained 
after 229 iterations to achieve a one part in 10^ change in the cost function in an iteration 
(we will use this same convergence criterion throughout), had Cmass = 0.13, Cckm = 4.65, 
Cscnc = 3 X 10~^, and Cparameter = 38.9, giving a total cost function C = 43.7. The values of 
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the parameters giving this fit are as follows, 



= 71.27, Qr, = 71.27 , 
g; = 0.811, 4 = 0.0201, g; = 0.00831 , 

= 0.00715, g^ = 0.00112, g^ = 0.000371 , 
5iR = 0.00269, 52R = 0.0340, S^r = -0.0367 , 
Su = 0.00074, S21 = 0.0027, S31 = -0.0034 , 
e = 150.8 degrees , 



I 0.1612 0.0477 0.0268 \ 

0.0144 0.00024 -0.0133 
V -0.0442 -0.0367 -0.1562/ 



/ 0.1660 0.1589 0.0189 \ 

0. 0.0180 -0.0190 
V -0.1398 -0.0189 -0.1841/ 



/ 0.1038 0. -0.0517\ 

[/?^] = 0. -0.00081 -0.0375 . (54) 
^-0.0366 -0.00011 0.0230 y 
We see that the largest value of the /3 asymmetry parameters is 0.184 in magnitude, so 
the first question we must address is whether this large asymmetry is needed to reproduce 
the large mixing S12 = 0.221 between the first and second families. To show that this is not 
the case, we exhibit the result of rerunning the fit, this time omitting the S13 and S23 terms 
from the cost function. The result, attained after 137 iterations, has Cmass = 0.04 (that is, 
the fitted mass values are right on their targets) and S12 = 0.221, so that the Cabibbo mixing 
is also right on target, but the largest of the /3 asymmetry parameters has a magnitude of 
0.01, a factor of 18 smaller than in the fit of Eq. (54). The values for the unconstrained third 
family mixings obtained this way are S13 = 0.00021, S23 = 0.00072, much smaller than in 
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the fit of Eq. (54). So we conclude that the large P asymmetry values of Eq. (54) are needed 
to get correct fits to the third family mixings; the correct value of S12 by itself is obtained 
with /3 values much smaller in magnitude than S12, in agreement with our observation in 
Sec. VII that S12 is of zeroth order in the asymmetries. 

As a second experiment, which gives further insight into why the model requires large 
asymmetries to fit the third family mixings, we rerun the fit replacing the targets for both 
Si3 and S23 by their geometric mean ~ 0.011, with a standard deviation of 0.0015. We find 
now convergence in 216 iterations, with Cmass ==0.02 (that is, again the fitted mass values 
are right on their targets), and values for the CKM mixings of S12 = 0.221, S13 = 0.0117, 
S23 = 0.0101. For the other components of the cost function we find Cckm = 0.86, Cgcnc — 
0.9 X 10~^, and Cparameter — 3.5, for a total of C = 4.4. As suggested by the small value 
of C'parameter, the largest of the /3 asymmetry parameters now has a magnitude of 0.028, a 
factor of 6.6 smaller than in the fit of Eq. (54). We conclude from this fit that what requires 
the large asymmetries in Eq. (54) is splitting S23 and S13 from a common mean value. 

This conclusion can be understood from a simple analytic model, in which the corrections 
of Appendix C to U[ and U{i are neglected. Referring to Eq. (42a), let us write Vckm to 
first order accuracy as 

/ 1 ^12\ 

Vckm = , (55a) 

SO that the fitted S12 is given by S12 = |t'i2|- Then from the approximation of Eq. (42b) for 
f^CKM, together with the CP invariance condition [see Eq. (38b) and the discussion following 
Eq. (39c)] (7(3 = (723, fi^d that 

Sl3 = \S3 - dsl/S, S23 |S3 + dsl/S , 

S3 = - (7^3, ^3 = -^120-^3 ■ (55b) 

Thus, the spread of S13 and S23 from their geometric mean is governed by ^3, in which the 
quantity (723, which is a linear combination of the P asymmetries, is suppressed in magnitude 
by a factor of 1^121 = S12 = 0.221. This is why large P asymmetries are needed to fit the 
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experimental data, whereas much smaller asymmetries suffice when the observed S13 and 
S23 are replaced in the fitting program by their geometric mean. For example, in the fit of 
Eq. (54), the magnitude of ds is 0.0445, which corresponds to a value (T23 = 0.0445/0.221 = 
0.20, similar in size to the maximum /3 asymmetries found in the fits. Thus, the six Higgs 
model interprets the large difference in magnitude between the observed S13 and S23 as 
indicating asymmetries in the Yukawa couplings substantially larger than one might naively 
infer from the magnitude of S23- The possible relevance of this observation to the extension 
of our model to neutrino mixing will be discussed in Sec. X. 

We next address issues of fine tuning and naturalness in the six Higgs model. In the fit 
of Eq. (54), the absolute values of the matrix elements of the matrices U'l''^ and take 
the values 

/ 0.974 0.224 0.055 \ 
[\U2\] = 0.224 0.974 0.034 , (56a) 
^0.046 0.045 l.OOOy 



\ut\] 



( 1.000 0.00010 0.066 \ 

0.00010 1.000 0.067 

V 0.066 0.067 1.000 y 



(56b) 



( 1.000 
0.00010 

V 0.033 



0.00010 0.033 \ 
1.000 0.036 

0.036 1.000 j 



(56c) 



We see that the mixing S12 of the first two families arises nearly entirely from [/^, while the 
2x2 submatrices of and C/^^, that mix the first two famihes (and that are equal to good 
accuracy) are nearly the unit matrix, which is what allows the strangeness changing neutral 
current constraint to be satisfied. To estimate the amount of fine tuning involved in this, we 
note that l-Dniggsl of Eq. (47b) is quadratic in the matrix element \U'ivi\ — I^^rLI- Hence if the 
entry for |i^Higgs|/ (2.6 x 10~^"^) in Table III were scaled up from 0.016 to unity, corresponding 
to the strangeness changing neutral current constraint being just barely satisfied, the off 
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diagonal matrix elements |C/fi2| — |t^fli2l i^i Eqs. (56b, c) would be scaled up from 0.00010 
to O.OOOIO/O.OI62 ~ 0.00079. Taking as a "generic" off diagonal matrix element the average 
value ~ 0.05 of the 13 and 23 matrix elements of Eqs. (56b, c), we estimate that fine tuning 
in the mixing matrices, of order a factor of 0.05/0.00079 ~ 63, is involved in satisfying the 
strangeness changing neutral current constraint, for an assumed Higgs mass in the fit [see 
the second hne in Eq. (46b)] of M^j"^ = (4.5 x 0.3)^0^ ~ 83 GeV. For a Higgs mass of 330 
GeV the fine tuning would be correspondingly reduced to a factor of roughly 16, and for a 
Higgs mass of 800 GeV the fine tuning factor would be roughly 6. 

Given that there is some fine tuning involved in obeying the strangeness changing neutral 
current constraint, one can ask whether it is natural or unnatural to the experimental data. 
If the fine tuning is not natural to the data being fit, one would expect the fits to the masses 
and CKM parameters to improve, or the convergence to a fit to become faster, when the 
cost function term Cgcnc is omitted from the total cost function. Performing this experiment, 
we find that without Cgcnc; a comparably good fit is obtained (Cmass = 0.64, Cckm = 4.1) 
as with the cost function term Cgcnc included, but 600 iterations, as opposed to 229, are 
required for comparable convergence. In other words, the strangeness changing neutral 
current constraint appears to guide the search to a region of parameter space that gives a 
good fit; we interpret this as an indication that the fine tuning involved in satisfying this 
constraint is in fact natural to the data. 

One other place where there is fine tuning in the fits is in the first family masses, since 
these are naturally zero only in the absence of Yukawa coupling asymmetries. In principle, 
if the first family cost function terms are omitted from Cmass; one might expect first family 
masses as large as 0.2 (the value of the maximum asymmetry parameters) times the corre- 
sponding third family masses, which would give ~ 35, ~ 0.9, Me ~ 0.4. However, 
performing the experiment of omitting first family mass constraints from the fit, we find 
first family masses M„ ~ 0.9, ~ 0.23 , Mg ~ 0.07, that is, the first family masses are still 
smaller (or in the case of M^, equal to) the second family masses. We interpret this as an 
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indication that small first family masses are in fact natural to the remaining experimental 
data when first family masses are excluded, in the framework of the six Higgs model. 

We conclude this section by giving some comparative fits in the three Higgs model. Using 
the same cost function parameters and convergence criterion as in the six Higgs case, we 
get the three Higgs model fit shown in Table IV, which required 864 iterations. The mass 
fit is generally good, except for the low value Ms = 0.037 (corresponding to Cmass — 7.3), 
while the CKM parameters are close to their targets (Cckm = 0.5). When the strangeness 
changing neutral current constraint is omitted in this case, we find faster convergence (398 
iterations) and a better fit, with Ms = 0.151 (corresponding to Cmass = 0.7) and with the 
CKM parameters right on target (Cckm = 0.1). This behavior contrasts sharply with what 
we saw in the six Higgs model fits, and we interpret it as indicating that the strangeness 
changing neutral current constraint is not natural to the data as interpreted in the three 
Higgs model. 

X. EXPERIMENTAL ISSUES, NEUTRINO MIXING, COUPLING CONSTANT 
UNIFICATION, AND DIRECTIONS FOR FUTURE WORK 

Of the two models that we have developed in the previous sections, we find the six Higgs 
model the more interesting as a candidate for an extension of the standard model into the 
energy region that will become accessible in the next decade. As compared with the three 
Higgs model, the six Higgs model has fewer parameters, gives better overall fits to the data, 
and gives some indication that the strangeness changing neutral current constraint is natural 
to the data. It also violates CP spontaneously in an interesting way that is correlated with 
the generation of second family masses for the u, d, e families. 

The prime experimental signature of the six Higgs model is the spectrum of Higgs states 
tabulated in Table II. If the potential V2 that couples the 4> to rj Higgs overall phases is in 
fact small, then the lightest Higgs states should be the pseudo-Goldstone bosons. However, 
because of the | power scahng law of Eq. (25c), they need not be so light as to conflict 
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with current Higgs mass limits. For example, if the Higgs masses Mniggs that enter into the 
strangeness changing neutral current constraint are of order 330 GeV, and V2/V^,r? ~ 0.1, 
which is in the weak coupling regime, then the pseudo-Goldstone boson masses are expected 



Although we have included the possibility of a right handed neutrino, and of Dirac 
neutrino masses and mixing analogous to CKM mixing, in our Lagrangian, we have not 
attempted a detailed study of the neutrino sector because the experimental picture there is 
still incomplete. However, let us briefly address the recent report by the Super-Kamiokande 
Collaboration [16] of evidence for atmospheric neutrino oscillations, suggesting large mixing 
(of order unity) of second and third family neutrinos. This is clearly a different pattern than 
is seen for the charged fermion mixings, where, for example, in the fit of Eq. (54) the ^ — r 
mixing matrix elements of t/£ are smaller than 0.01 in magnitude. Large i/^ — i^r mixing can 
be accommodated in our model, nonetheless, by assuming that the Yukawa coupling ratio 
g^/gjp, which we have taken to be small for f — u,d, e, is close to unity for f — v. Together 
with Vtr}/^<t> — 1, this implies that Ry of Eq. (39b) is close to unity in magnitude (although 
it can have a nonzero phase). Referring to Eqs. (38a, b), we see that this imphes that the 
neutrino mass matrix is now nearly degenerate in the two dimensional subspace spanning 
the second and third families, and so small asymmetries, or asymmetries nearly equal in 
magnitude, then imply nearly maximal mixing. To show this explicitly, let us apply the 
analysis of Appendix B to the mass matrix 



to be of order (0.1)2 330 GeV ~ 104 GeV, above current Higgs mass limits. 




(57a) 



Then for Ml — mw) , we have from Eq. (B2b), 




(57b) 



and so Eq. (B4b) gives for the mixing angle 




(57c) 
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Thus there is maximal mixing whenever 

^|i?V3*2 + (J23|»|i?1'-l + ^(|<J23r-k32r) ■ (57d) 

If lit!*^! is close to unity, this inequality can be satisfied either (i) if (723 and (732 are both 
small, or (ii) if the magnitudes of (J23 and (T32 are not small, but are approximately equal. 
In Sec. IX we saw that to reproduce the observed CKM parameters S23 and S13, we needed 
sizable asymmetries (of order 0.2), which if also present in the neutrino sector /3's would 
allow near maximal mixing of the second and third family neutrinos by case (ii) even when 
the ratio R,y is only approximately unity in magnitude. Thus large — mixing is easy to 
achieve in the six Higgs model. Less natural is near degeneracy of the masses of Ve and i/^, 
as appears to be needed for both the MSW and the vacuum oscillation interpretations of 
the solar neutrino data, since the first family masses are zero in our model in the absence of 
Yukawa asymmetries. Such a degeneracy would have to be the result of sizable asymmetries 
together with substantial fine tuning in the neutrino mass matrix, either to raise the mass 
to close to the Ufj, mass in case (i), or to lower the mass to close to zero in case (ii) (as, 
for example, is done in the model of Barger, Pakvasa, Weiler, and Whisnant [17]; see also 
Baltz, Goldhaber, and Goldhaber [17]). In either case, there will almost certainly be large 
mixing of Pe with i/^; so on this (very preliminary) interpretation, our model would favor 
the large angle as opposed to the small angle MSW solution. 

Let us next address the issue of coupling constant unification in the six Higgs model. 
Because we do not alter the fermion representation content of the standard model, the 
usual running coupling analysis applies. As noted by Langacker [18], the standard model 
with ~ 7 (by current data [19], 7.66) Higgs boson doublets gives one loop coupling constant 
unification with a unification energy of order 5 x 10^^ GeV. Even with only 6 Higgs doublets, 
the magnitude of two loop radiative corrections [20] is sufficient to make coupling constant 
unification a possibility. Of course, because the unification energy is lower than in the 
customary scenario, a mechanism is needed to suppress proton decay, such as is present 
in the SU{15) family [21] of grand unification models. Clearly, definitive statements here 
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will depend on the nature of the high energy theory for which the six Higgs model is a low 
energy effective theory; the point we wish to stress, though, is that the six Higgs model 
may be a candidate for coupling constant unification without the assumption of low energy 
supersymmetry. Whether such a candidate is needed, of course, will depend on the outcome 
of supersymmetry searches over the next decade. 

There are a number of obvious directions for further work on the models we have de- 
veloped in this paper. Entirely within the low energy effective action framework, one can 
address the issue of one loop radiative corrections to the mass and mixing matrix analysis 
given here. This will involve the parameters determining the Higgs masses in an integral way, 
and if the six Higgs model is to be viable, the one loop corrections should improve, rather 
than make worse, the comparisons with experiment and the consistency tests discussed in 
Sec. IX. Another issue that can be addressed within the low energy framework is the mag- 
nitude of electroweak baryogenesis in the six Higgs model, and cosmological implications of 
this model more generally. At a deeper level, there is the issue of finding a grand unified 
model, composite model, or hybrid model comprising elements of both, which is a natural 
high energy physics source for the low energy effective action physics described by the six 
Higgs model. Such a high energy model must, through its representation content and in- 
stanton physics, justify the discrete chiral transformation rules assumed in Eqs. (la-c), and 
it is also the place where one must seek explanations for the "vertical" hierarchy of Yukawa 
coupling strengths, and the pattern of Yukawa coupling asymmetries, that is needed for our 
fits. 
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Added note. After this paper was posted to the Los Alamos e-print archive, two earlier 
papers that use families of Higgs scalars (although without the ingredient of Zq discrete chiral 
symmetry analyzed here) were brought to my attention. The paper of Derman and Jones 
[22] studies a two family, two Higgs doublet model with an S2 permutation symmetry, and 
is probably the earliest paper to extend the idea of family symmetries to the Higgs sector; 
the paper of Derman [23] extends this to three famihes of fermions and Higgs doublets with 
an Ss permutation symmetry. 

APPENDIX A: NUMERICAL MINIMIZATION OF THE HIGGS POTENTIAL 

Because the Higgs potential of Eqs. (7a-c) is comphcated, even with the simplifying 
assumptions of CP invariance and cyclic permutation symmetry, we have supplemented 
our analytic studies of the Higgs extrema with numerical studies, performed by using the 
conjugate gradient method to minimize the Higgs potential. Since it is easy to analytically 
compute the first derivatives (the gradients) of the Higgs potential, it is advantageous to 
use the conjugate gradient method in a form where both the function to be minimized and 
its derivatives are externally supplied; this gives a faster routine and there is some built in 
redundancy that serves as a check, since the same information is in effect furnished twice, 
once through the computation of the function and a second time through the independent 
computation of its derivatives. We have used the minimization program "frprmn" of Press 
et al. [12], with the following modification. Press et al. base the convergence criterion in 
their program on computing the change in the function value over one iteration, but this 
results in significant truncation error inaccuracies for the minimizing values of the arguments 
(the Higgs fields) when the function is large in magnitude but very fiat at its minimum. 
Since the gradients are explicitly known, and since at the minimum the gradients must 
all vanish, much better accuracy for the minimizing Higgs fields is obtained by making the 
convergence criterion depend on the maximum gradient. With this modification to "frprmn" , 
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one can verify vanishing of the gradients to double precision accuracy at the Higgs potential 
minimum. 

To obtain the formulas for the gradients, as a function of general 4>, r], we substitute 
(f)^(f)-\-5(f),r]^r]-\-5r] into Eqs. (7a-c), and retain the first order variations, which can be 
brought to the convenient form 

^-^Higgs potential — 

Re^(C„^50; + CM . (Ala) 

n=l 

We assume both CP invariance and cyclic permutation symmetry; using the latter we get 
formulas for Cp^ by cychc permutation of the arguments of Cf'''. Changing notation for 
the coefficients from C^mn to Cg-mn, to avoid notational ambiguities when explicit numerical 
values are assigned for m, we obtain the following explicit expressions for Cf'^. 

Cf = 4A^(0t01 - vl)(j)l - 2(/ii0 + /X20)(0202 + 03^3)01 

-Q;^ (2020^03 + 03</'2 + 0203) + IVl 

3 

+ [2Ci 

m=l 

+C3.im4>2r)mVm-l + C3.3m4>3Vm-lVm 
+C4;ml(l>3VmVm+l + C'4;m202??m^m+1 
+C5;lmr]2r}*m<Pm-l + Cr,.rn2m^m+l<Pm 
+C6;ml^3^^0m+1 + Ce-SmVaVm-l^t^m 

+C8;lm02^^0m-1 + C'8;3m03?7m0m-l + C'8;m20m0^+1^2 

■,m2VmVm+lV2 
+Cu;lmr]2(p*mr)m-1 + C'll;ml0m Wl^3 

+Cl2;3m?73?7m0m-1 + C'i2;m2?7m0m+1?72] , (Alb) 
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-arj{2r]2vlV3 + vlv*2 + vWa) + 70i 

3 

+ Yl [2Cl;ml0i^0m?7l 

m=l 

+C2-ml<Pl<P*mVm + C'2;lm010^^7m 
+C3;ml%0m^m+1 + Ci;m2r]2<t>m<Kn+l 
+CA-lmV2<Pm4>m-l + C'4;3m'730m0m-1 
+C5;ml030m?7m+1 + C^.3rn(t>3(Kn-lVm 
+C6-lm<p2(p*mVm-l + C'6;m2020^+ir7m 
+C7;m20m0i^+102 
0m+103 

+C9;lm?720m^m-1 + C'9;3m?730m?7m-1 + C'9;m2?7m?7m+1^2 
+ClO;lm^2»7;^0m-l + C'lO;3m%'7m0m-l + C'lO;ml»7m Wl03 
+Cll;3m030m?7m-1 + C'll;m20m?7m+102 

+Cl2;lm02?7j^^m-1 + C'l2;ml'7j^0m+103] ■ 

(Ale) 



APPENDIX B: BI-UNITARY DIAGONALIZATION OF A 2 x 2 MATRIX 

We give here the method for constructing the matrices Vl and that obey Eq. (41a) 
of the text, suppressing the flavor index / throughout. Let m by the 2x2 complex matrix 
defined by 

(<7li cri2\ 
(Bl) 
(721 (722/ 

We begin by forming the self-adjoint matrices Ml = mm^ and Mr = m^m, which we write 
in the form 

Ml = {"^"^ \ , Mn^l '], (B2a) 
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with 



I |2 I I |2 



Bl = |cr2lP + |cr22|^ , 

Zl = crl^cr2i + oX-^P^i-i > (B2b) 



and 



A-R — \(T\\f + |(72l|^ , 

Br = |cri2|^ + |cr22p , 

^i? = Cri2'^ll + C^22<^21 ■ (B2c) 

The quantities just defined are not independent, since it is easy to verify that 

Al + Bl^Ar + Br , 
^{Al - Bl)' + \zl\' = \{Ar - BrY + \zr\' , 

\zl?<AlBl, \zr\'<ArBr . (B2d) 

The desired bi-unitary matrices will be the Vl for which VlM^vI is diagonal, and the Vr 
for which VrMrV^ is diagonal, with eigenvalues ordered in magnitude. 
Thus, defining the self-adjoint matrix M by 

/A z*\ 

M = \ , \z\'' < AB , (B3a) 

\z B) 

it suffices to find the diagonalizing unitary transformation V that yields 

\ 

VMV^ = , < \k2\ ; (B3b) 

V |«2|V 

then all that we have to do is to apply this construction twice, first to Ml and then to Mr. 
Let us write M in Pauli matrix form as 

M^]^{A + B)+v-T, v^[zR,zi,]^{A-B)j , (B3c) 
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with zrj the real and imaginary parts of z. Representing the diagonahzing V in Pauh 
matrix form as 

V — exp(i©ri • t) = cos © + m • r sin © , (B4a) 

and letting z = (0, 0, 1) be the unit vector in the third axis direction, a simple calculation 
shows that we satisfy Eq. (B3b) by taking 

sin 2© = 



\Z X v\ 



\v\ 



[l(A-Br + \z\t- 



cos 2© = —— = ^ ' 



\v\ 



^1 1 / -2UI 
© = -tan"^ 



and that this V gives 



2 V^-^y 

\z X v\ \z\ 



VMV^ = ^(A + - \v\t^ . (B4c) 

Thus we see that the squared eigenvalues are 

= + = + + , (B5a) 

which are correctly ordered; the smaller squared eigenvalue is guaranteed to be nonnegative 
by virtue of the fact that the product of the squared eigenvalues is 

^-^{A^Bf -\v\^ = AB -\z\^>^ . (B5b) 

When \z\ =0, the above formulas are indeterminate; we then get the correct eigenvalue 
ordering by taking sin 2© = and cos 2© = ±1, with the + sign holding for A < B and the 
— sign holding for A > B. Referring back to the identities of Eq. (B2d), we see that they 
imply that \vl\ = \vr\, and thus the eigenvalues are the same for Ml and Mr, as expected. 
Substituting the expression for n in Eq. (B4b) back into Eq. (B4a), we get the further useful 
expression 

(cos© — frsin©" 
1^1 I . (B6) 

1^ sin © cos © 
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APPENDIX C: IMPROVED FORMULAS FOR U[ AND Uj^ 



f ft 

In our numerical work, we used an improved approximation to Ul and C/^ obtained by 
adding to Eqs. (41b, c) the respective corrections AU[ and AUj^ , given by 













1 / /* 



\ '132^ L21 9^*2 '132^ L22 



lrf*'r,f T// 
9^ 



(Cla) 



and 



/O 


Vo 








l^R12'^2 ^23 \ 

lT//t u;f*nf 
9 '^R22'^2 '123 



iff* n 
■9'^2^23 ^ 



(Clb) 



Here V/^ are the matrices defined in Eq. (41a) and computed in Appendix B, is the 
eigenvalue defined in Eq. (41a), given explicitly by 



'^2 — ^L2\{^il^R\2 + "^12^22) + ^22 ("^21^^12 + '^22^^22 



./ T//t 



/ (^S T//t 



./ T//t 



and the quantities r^ls) ''732 ^-^^ defined by 



^723 ~ ^/2l'^13 + ^/22'^23 ) 



(C2a) 



^32 — "^31^12 + '^32^22 



(C2b) 



These corrections make the formulas for U[ and Uj^ accurate to first order when (|«;^|/3)^, 
rather than 1^21/3, is regarded as a first order small quantity. They have only a small effect 
on the fits of Sec. IX (because for charged fermions the second to third generation mass 
ratios are small), but are useful in performing accurate numerical checks that U{^M'fUQ is 
diagonal. 
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TABLES 



TABLE I. Higgs eigenmodes, masses, and fermion couplings for the 3 Higgs doublet model in 
the cyclic symmetry limit 



mode 
designation 


charge 


mass 

squared 


fermion family 
couplings 


^R '^7 
^R 


±1 
±1 








3{fj.2 + a) 02 
3(/X2 + a)n'^ 

(4A + 2ni + 2/X2 + la)n^ 
|aJ^2 


1st 

2nd 

3rd 
1st and 2nd 
1st and 2nd 


TABLE IL Higgs eigenmodes, masses, and fermion couplings for the 6 Higj 
the cyclic symmetry limit, assuming weak coupling of (/> to r/ 


IS doublet model in 


mode 
designation 


charge 


mass 
squared 


fermion family 
couplings 


,(3) 
^<t>R 
,(+) g(-) 

^^R 


±1 
±1 










1st 

2nd 

3rd 
1st and 2nd 
1st and 2nd 


r(2) 
Or, 

^r,R 

1+) A-) 

J-) 

^TjR ^^vl 


±1 
±1 






3(/X2„ + Q;^)n2 
4A^v2 

(4A^ + + 2//2^ + larj)n'^ 
1^2 


3rd 
1st 
2nd 
1st and 3rd 

1st and 3rd 


r(3) 
"PG 
.(3) 
^PG 


±1 






2nd and 3rd 
2nd and 3rd 
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TABLE III. Six Higgs model fit to experimental data 



quantity 


target value 


fitted value 




247. 


247. 




0. 


0.001 


Mu 


0.005 


0.005 


M, 


1.30 


1.28 


Ml 


173. 


173. 


Md 


0.010 


0.011 


Ms 


0.200 


0.219 


Mb 


4.30 


4.29 


Me 


0.00051 


0.00051 




0.1057 


0.1057 


Mr 


1.777 


1.777 


Pffiggsl 

2.6x10-1* 


0. 


0.016 


Sl2 


0.221 


0.221 


Sl3 


0.0035 


0.0041 


S23 


0.041 


0.035 


|sin (5i3| 


0.60 


0.44 



TABLE IV. Three Higgs model fit to experimental data 



quantity 


target value 


fitted value 


V = 62O0 


247. 


247. 


M„ 


0.005 


0.005 


Me 


1.30 


1.28 


Mt 


173. 


173. 


Md 


0.010 


0.011 


Ms 


0.200 


0.037 


Mb 


4.30 


4.31 


Me 


0.00051 


0.00051 


M^ 


0.1057 


0.1057 


Mr 


1.777 


1.777 


I^Higgsl 

2.6x10-1-1 


0. 


0.001 


Sl2 


0.221 


0.221 


Sl3 


0.0035 


0.0037 


S23 


0.041 


0.039 


|sin Sisl 


0.60 


0.55 
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